AD-A229  138 


OTIC  FILE  COPY 


NASA  Contractor  Report  187^41 
ICASE  Report  No.  90-68 


ICASE 


REGULARIZED  CHAPMAN-ENSKOG  EXPANSION 
FOR  SCALAR  CONSERVATION  LAWS 


Steven  Schochet 
Eitan  Tadmor 


Contract  No.  NAS  1-18605 
October  1990 


DTIC 

ELECTE 


DECO  3 1990 


B 


Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  Virginia  23665-5225 

Operated  by  the  Universities  Space  Research  Association 


fVI/\S/\ 

National  Apronaiilirs  and 
Space  Administration 

Langley  Research  Center 

Hampton,  Vifginia  23665  5?25 

""dTsTRIBUTION  STATEMarr  IT 

Ap  proved  for  public  rdbaMJ 
r.stdbation  Ualimltedi  . 


REGULARIZED  CHAPMAN-ENSKOG  EXPANSION 
FOR  SCALAR  CONSERVATION  LAWS 


Steven  Schochet  and  Eitan  Tadmor^ 

School  of  Mathematical  Sciences 
Raymond  and  Beverly  Sackler  Faculty  of  Exact  Sciences 
Tel  Aviv  University 
Tel  Aviv  69978  Israel 


^  ABSTRACT 

Rosenau  [Phys.  Rev.  A,  40  (1989),  pp.  7193-6]  has  recently  proposed  a  regularized 
version  of  the  Chapman-Enskog  expansion  of  hydrodynamics.  This  regularized  expansion 
resembles  the  usual  Navier-Stokes  viscosity  terms  at  law  wave- numbers,  but  unlike  the  latter, 

it  has  the  advantage  of  being  a  bounded  macroscopic  approximation  to  the  linearized  collision 

r'. 

operator.  7  f ^  i  I  ,  ?  ‘ 

Q. _ In  this  paper^fC'  study^e  behavior  of  Rosenau  regularization  of  me  Chapman-Enskog 


expansion  (R-C-E)  in  the  context  of  scalar  conservation  laws.  We-afeow-that  this  R-C-E 
model  retains  the  essential  properties  of  the  usual  viscosity  approximation,  e.g.,  existence 
of  travelling  waves,  monotonicity,  upper-Lipschitz  continuity  etc.,  and  at  the  same  time,  it 
sharpens  the  standard  viscous  shock  layers.  We  prove  that  the  regularized  R-C-E  approx¬ 
imation  converges  to  the  underlying  inviscid  entropy  solution  as  its  mean-free-path  ^10, 
and  we  estimate  the  convergence  rate.  A  ,  \  /  \ 


on 
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1.  Introduction 


Rosenau  [R]  has  recently  proposed  the  scalar  equation 

ct  \  [ 

Uf  +  /(u)x  =  — 

(1.1) 


T  V 


=  £ 


1  + 

1 


u{k) 


-u{k) 


Ll+mVP 

as  a  model  for  his  regularized  version  of  the  Chapman- Enskog  expansion  for  hydrodynamics. 
The  operator  on  the  right  side  looks  like  the  usual  viscosity  term  su^x  at  low  wave-numbers 
k,  while  for  higher  wave  numbers  it  is  intended  to  model  a  bounded  approximation  of  a 
linearized  collision  operator,  thereby  avoiding  the  artificial  instabilities  that  occur  when  the 
Chapman-Enskog  expansion  for  such  an  operator  is  truncated  after  a  finite  number  of  terms 
[R]. 

In  this  ajticle  we  will  compare  the  behavior  of  solutions  of  Rosenau  regularization  for  the 
Chapman-Enskog  expansion,  abbreviated  hereafter  as  the  R-C-E  equation,  (1.1),  with  those 
of  the  viscous  conservation  law 


(1.2) 


“1  +  /(“)•  = 


towards  which  (1.1)  tends  as  m  0,  and  with  the  conservation  law  with  absorption 
(1.3)  ut  +  f{u)x^ 


u 


£m‘ 


Since  the  right  side  of  (1.1)  tends  to  that  of  (1.3)  in  the  limit  of  large  fc,  it  is  not  surprising 
that  the  smoothness  properties  of  solutions  of  the  former  resemble  those  of  the  latter.  In 
particular,  the  R-C-E  equation  does  not  smooth  out  initial  discontinuities,  but  as  shown 
in  §2,  it  does  preserve  the  smoothness  of  initially  smooth  small  initial  data.  On  the  other 
hand,  the  right  side  of  (1.1)  also  resembles  that  of  (1.2)  in  that  both  are  second  derivatives. 
Consequently,  it  is  plausible  that  the  regularized  R-C-E  equation  (1.1),  like  the  ordinary 
viscosity  equation  (1.2),  should  have  travelling  wave  solutions  connecting  shock  states  of  the 
underlying  conservation  law 


(1.4) 


“i  +  /(“)•  =  0 


In  §3  we  show  that  when  f"  >  0  such  solutions  exist  iff  m  is  sufficiently  small. 

At  the  same  time,  solutions  of  the  R-C-E  equation  (1.1)  also  resemble  those  of  the 
inviscid  conservation  law  (1.4)  in  that  both  admit  unique  entropy  solutions  which  share 
similar  properties.  In  §4  we  show  that  the  R-C-E  solution  operator  associated  with  (1.1), 
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1 _ 


'.G3 


□  □ 


like  the  entropy  solution  operator  of  (1.4),  is  L^-contractive,  monotone,  and  BV-bounded. 
Furthermore,  the  R-C-E  solution  of  (1.1)  tends  to  the  inviscid  entropy  solution  of  (1.4)  as 
the  ’mean-free-path’  £  |  0.  Finally,  if  f"  >  0,  the  R-C-E  entropy  solution  of  (1.1)  is  also 
upper-Lipschitz  continuous,  in  agreement  with  Oleinik’s  E-condition  which  characterizes  the 
entropy  solution  of  (1.4),  and  in  §5  we  estimate  the  convergence  rate  of  the  former  to  the 
latter  as  £  J.  0. 

2.  Smoothness 

It  is  well  known  that  solutions  of  (1.2)  are  smooth  for  t  >  0;  i.e.  initial  discontinuities  are 
smoothed  out  at  positive  times.  In  contrast,  by  looking  at  piece-wise  constant  initial  data 
or  at  the  linear  case  /(u)  =  u,  one  sees  that  initial  discontinuities  of  solutions  of  (1.3)  are 
merely  attenuated,  not  smoothed  out,  at  positive  times.  Since  the  damping  of  (1.1)  is  less 
than  that  of  (1.3),  it  is  clear  that  (1.1)  also  does  not  smooth  out  initial  discontinuities.  On 
the  other  hand,  if  the  (£-independent)  initial  data  for  (1.3)  is  smooth  then  it  will  remain  so 
provided  that  m  is  sufficiently  small  (see  below).  The  next  Theorem  tells  us  that  the  same 
holds  for  the  R-C-E  equation  (1.1). 

Theorem  2.1.  The  solution  of  the  R-C-E  equation  (1.1)  remains  as  smooth  as  its  initial 
data, 

(2.1)  u(x,0)  =  ■uo(a;), 
provided  the  initial  data  Uq  are  sufficiently  small  so  that 

(2.2)  2{m||uo||i,oo||/''(iio)||Loo}^/^ -Frn^£||/"(uo)||z."o||uo||i«>  <  1. 

Remark.  Since  we  can  ensure  that  (2.2)  will  be  satisfied  for  any  fixed  initial  data  by 
making  m  sufficiently  small,  Theorem  2.1  can  also  be  viewed  as  showing  how  the  smoothness 
properties  of  the  R-C-E  equation  (1.1)  approach  those  of  the  viscosity  equation  (1.2)  as 
771  ^  0. 

Proof.  We  will  show  formally  that  (2.2)  implies  a  bound  on  the  L°°  norms  of  u  and  Ux. 
Estimates  for  higher  derivatives  then  follow  in  standard  fashion;  see  [M].  Furthermore,  this 
fact  ensures  that  the  formal  estimates  can  be  justified  either  by  smoothing  the  initial  data 
or  by  applying  a  further  vanishing  viscosity  regularization. 
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The  first  step  towards  obtaining  the  desired  bounds  is  to  note  that  the  right  side  of  (1.1) 


can  be  written  as 


-£P 
1  + 


=  -  lu—  - ; — u(k) 

m^e\  Ll+m^e^A:^ 


'{u  Qme  *  ) 


where  *  denotes  convolution,  and 


satisfies 


Q.(x)  =  le-l'l" 


=  1  • 


To  obteiin  a  uniform  bound  on  u,  multiply  (1.1)  by  and  integrate  over  x;  since 

|u^“^|ii/('u)j;  is  an  exact  derivative  its  integral  vanishes,  while  the  contribution  of  the  right 
hand-side  (2.3)  is  nonpositive,  for  by  (2.5), 

f  -  Qrne  *  u}dx  <  — ||u||^p  -  ||u||^;^  ||  } 

J  m*£  m^e 


m*£ 

Dividing  the  remaining  inequality  by  (p  —  l)||u||£,p  and  integrating  over  t  from  0  to  T,  we 
obtain 

(2.7)  llu(r)|Up  <  lh:o|U., 

and  the  boundedness  of  ||■u(^)llx,oo  follows  by  letting  p  t  oo. 

In  order  to  estimate  in  similar  fashion  the  L°°  norm  of  u*,  we  differentiate  (1.1),  obtaining 

(2.8)  u^t  +  f{y-)xx  =  — -  Q'm,  *  , 

TTh  c 

and  as  before,  we  multiply  (2.8)  by  \ux\^~'^Ux  and  integrate  over  x.  Integrating  by  parts 
where  necessary  in  the  term  containing  /,  and  noting  that  by  (2.4) 

(2.9)  ll<3'„.IU'  = 

we  obtain  after  factoring  out  (p  —  l)||uj;||£,p. 


(2.10) 


CLl  TTh  c 
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Next  we  denote 


(2.11)  Y{T)  =  m^£ll/>o)llL»lK(T)lk». 

Applying  Gronwall’s  lemma  to  (2.10)  and  letting  p  t  oo,  we  obtain  that  Y  =  supyy'(!r) 
does  not  exceed 

Y{T)<Y  <  + 

(2.12) 

<  Y(0)  +  ^\\uoU~\\f"{uo)\U., 

as  long  as  y  <1.  Estimate  (2.12)  is  a  quadratic  inequality  for  Y  for  which  the  roots  of  the 
corresponding  equation  are 

(2.13)  K  =  i{l  +  r(0)  ±  s/(l  -  Y(0W  -  4m||u„|U.n/"(u„)|U-}  . 

Our  assumption  (2.2)  tells  us  that  the  expression  under  the  square  root  on  the  right  is 
positive.  Since  y’(O)  is  bounded  by  the  smaller  root  in  (2.13),  it  follows  from  (2.12)  plus 
the  continuity  of  Y(T)  that  Y(T)  remains  bounded  by  this  root.  This  in  turn  confirms  that 
indeed  Y  <  1,  and  the  uniform  bound  of  j|ua;(7’)||£oo  follows.  □ 

Remark.  Arguing  along  the  above  lines  for  the  conservation  law  with  absorption  (1.3), 
one  arrives  at  the  inequality,  analogous  to  (2.12),  Y{T)  <  This  shows 

that  if  y'(O)  <  1  then  Y{T)  remains  <  1.  Consequently,  if  y’(O)  <  1,  then  Y{T)  e'lj  hence 
[[Ui(T)|(£,<»,  satisfy  a  maximum  principle  in  this  case. 

3.  Shock  Profiles 

Lax’s  generalized  entropy  conditions  [L]  for  ’’legitimate”  shock-wave  solutions  of  the 
conservation  law  (1.4)  can  be  interpreted  as  the  requirement  that  these  shocks  can  be  realized 
as  the  limit  of  travelling  wave  solutions  of  the  viscosity  equation  (1.2).  If  the  flux  function 
/  is  convex,  these  conditions  reduce  to  the  shock  inequalities  [L] 

(3.1)  /'(u_)>5>/'(u+), 

where  s  is  the  speed  of  the  shock  joining  on  the  left  to  on  the  right.  In  this  section 
we  show  the  analogous  result  for  the  (convex)  R-C-E  equation  (1.1):  It  admits  travelling 
wave  solutions  whose  limit  as  e  i  0  are  shock  wave  solutions  of  (1.4),  iff  (3.1)  holds  and  m 
is  sufficiently  small. 

Theorem  3.1.  Assume  f"  >  0.  Then  (3.1)  and  the  Rankine-Hugoniot  shock  condition 

(3.2)  //■(u+)  =  0,  ^(w)  = -3{u  -  u_} -I- {/(u)  - /(u_)}. 
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are  necessary  conditions  for  the  existence  of  a  travelling  wave  solution 

u{z  =  —j—),  lim,_±oou(z)  =  u±, 

for  (1.1). 

Conversely,  if  (3.1),  (S. 2)  hold,  then  a  sufficient  condition  on  m  for  the  existence  of  such  a 
travelling  wave  is 

(3.3)  4m^  sup  {—f"{u)H{u)}  <  1  , 

u+  <u<u_ 

and  a  necessary  condition  is 

(3.4)  <  1  . 


Here  u,  is  defined  by 

(3.5) 


f'{u,)  =  s  . 


Proof.  Define  z  —  and  let  '  denote  Using  (2.4)  we  find  that  a  solution  of  (1.1)  of 
the  form  u  =  u{z)  satisfies 

(3.6)  -su'  +  f{uy  =  {Qm  *  u}"  , 

where  the  convolution  is  now  taken  w.r.t.  the  variable  z.  The  condition  limx.+ioo'i^ 
implies  that  also  Qm  *  u  u±  as  z  ±oo,  so  there  exist  a  sequence  of  values  zf  tending 
to  ±oo  on  which  {Qm  *  'u)'  tends  to  zero.  Hence,  integrating  (3.6)  from  z~  to  z  and  letting 
ji  — ♦  oo  we  obtain 

(3.7)  H(u)  =  -3{u  -  U_}  +  {f(u)  -  /(Ti_)}  =  {Qm  *  u}'  . 

Now  letting  z  tend  to  +oo  along  the  sequence  we  find  from  (3.7)  that  i/(u+)  =  0,  i.e. 
(3.2)  holds. 

Noting  that  H"  =  /"  >  0,  we  see  that  H{u_)  =  0  =  H{u+)  implies  H'  =  f  —  s  <  0  at 
the  lessor  of  u±,  and  H"  >  0  at  the  greater  of  the  two.  Hence,  if  <  u_  then  (3.1)  holds, 
while  if  this  inequality  is  reversed  then  so  are  those  of  (3.1),  i.e.,  we  can  replace  (3.1)  by  the 
condition 

(3.8)  u_  >  . 
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Next,  we  apply  to  (3.7)  the  operator  1  —  (the  inverse  of  the  operator  of  convolution 

with  Qm),  to  obtain 

(3.9)  u'  =  {1  -  =  H{u)  -  m\H'{u)u"  +  H"{u){u'Y}  . 

dz^ 

We  note  that  since  all  nonzero  solutions  5  of  {1  —  =  0  are  unbounded  on  R,  the 

solution  of  (3.9)  with  bounded  u  and  u'  which  we  construct  below,  also  satisfies  (3.7).  To 
construct  such  a  solution  we  introduce  the  auxiliary  variable 

(3.10)  = 
which  enables  us  to  rewrite  (3.9)  as  the  2x2  system 

(3.11)  u'^v 

(3.12)  m^H'(u)v'  =  H(u)  —  v  —  H"{u)v^  . 

The  convexity  of  H(u)  together  with  the  Rankine-Hugoniot  condition  (3.2)  imply  that  the 
only  critical  points  of  system  (3.11-12)  are  (tt_,0)  and  (u^.,0). 

We  remark  that  the  linearization  of  (3.11-12)  near  the  critical  points  (u_,0)  and  (u+,0) 
shows  that  they  are  both  saddles,  so  that  topological  methods  (see  e.g.  [S])  cannot  be 
applied;  one  might  even  tempted  to  conclude  that  the  existence  of  a  trajectory  joining  these 
saddle  points  is  unlikely.  What  saves  the  day,  however,  is  the  fact  that  the  system  is  singular 
on  the  line  u  =  u,,  i.e.,  that  the  coefficient  H'{u)  on  the  left  of  (3.12)  vanishes  at  u,,  which 
by  (3.1)  lies  between  u_  and  u+. 

The  key  to  finding  a  trajectory  joining  the  two  critical  points  is  to  note  that  solutions  of 
(3.11-12)  can  cross  the  line  u  =  u,  only  at  points  (u.,v,)  which  make  the  right  side  of 

(3.12)  vanish:  Equation  (3.11)  implies  that  H'{u{z))  is  0{z  —  z,)  near  the  value  2*  for  which 
u{z,)  =  ti,,  and  hence  (3.12)  shows  that  |u|  — »  00  as  z  — »  z,,  unless  the  right  side  of  (3.12) 
tends  to  zero.  Also,  since  the  right  side  of  (3.12)  is  quadratic  in  v,  a  comparison  of  (3.12) 
with  the  equations  zv'  =  shows  that  in  fact  |ul  reaches  infinity  before  u  reaches  u,. 

In  order  to  obtain  a  trajectory  joining  u_  at  z  =  —00  to  at  z  =  -foo  it  is  therefore 
necessary  and  sufficient  to  find  trajectories  joining  (u-,0)  (respectively,  (u+,0))  at  z  =  —00 
(respectively,  z  =  -|-oo)  to  (u,,u,)  (with  the  same  value  of  u,  for  both  cases)  at  some  finite 
values  of  z;  we  can  always  arrange  for  the  two  values  of  z  to  coincide  because  the  system  is 
autonomous.  Since  trajectories  through  (ii*,t;*)  are  not  unique,  the  existence  of  our  desired 
trajectories,  which  when  put  together  join  u_  to  no  longer  seems  so  unlikely. 
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Now  the  right  side  of  (3.12)  is  a  quadratic  expression  in  v,  whose  roots  are 


(3.13) 


U±(ti) 


—  1  ±  ^1  + 


If  the  argument  under  the  square  root  is  negative  at  u  =  u*  then  clearly  no  such  u,  exists; 
this  gives  the  necessity  of  (3.4)  for  the  existence  of  travelling  wave  solutions.  We  now  turn 
to  discuss  the  sufficiency  of  condition  (3.3):  it  says  that  exist  for  all  u  between  u_  and 
u+,  and  we  want  to  show  that  when  this  happens,  then  the  trajectories  mentioned  above 
exist  iff  (3.8)  holds.  Namely,  that  these  trajectories  exist  if  we  replace  u_  and  by 


(3.14) 


u  =  maoc{u_,ii+}  and  u'^  =  min{'u_,u+}  , 


but  not  with  these  latter  two  interchanged. 

The  linearized  system  around  the  two  critical  points  has  the  form 


(3.15) 


Since  the  determinant  of  the  matrix  on  the  right  of  (3.15)  is  negative,  both  critical  points 
are  saddles,  as  claimed.  Now,  it  is  not  hard  to  calculate  directly  the  asymptotic  directions 
of  the  solutions  that  approach  each  critical  point  as  z  tends  to  ±oo,  as  these  are  simply 
the  eigenvectors  of  the  matrix  in  (3.15),  but  in  any  case  we  will  have  to  determine  from 
(2.11-12)  the  signs  of  u'  and  v'  in  various  regions,  and  this  information  suffices  to  determine 
which  region  the  various  asymptotic  directions  lie  in.  In  this  way  we  obtain  Diagram  1,  the 
phase-plane  diagram  for  the  case  when  m  satisfies  (3.3). 

Based  on  Diagram  1,  the  existence  of  a  travelling  wave  solution  is  argued  as  follows.  There 
is  a  Lrajectory  that  leaves  the  critical  point  (^“,0)  and  enters  region  I.  If  this  trajectory 
remains  in  region  I  until  u  reaches  the  value  it*,  then  by  the  above  analysis  it  reaches  the 
point  (it,,  u,);  in  this  case,  it'  =  u  as  well  as  it  are  monotonic  on  this  semi- trajectory.  The  only 
way  that  the  trajectory  can  leave  region  I  before  reaching  the  line  it  =  it,,  is  by  entering 
region  77;  but  u'  >  0  in  this  region,  so  “clearly”  the  trajectory  still  reaches  (it,,u,).  A 
similar  analysis  backwards  in  the  “time”  z  shows  that  there  is  a  semi- trajectory  from  (u,,u,) 
to  (ii+,0).  By  checking  the  other  trajectories  leaving  and  entering  each  critical  point  we  see 
that  no  trajectory  joins  u'*'  to  u~  or  either  point  to  itself. 

Although  the  above  argument  is  sound  provided  that  Diagram  1  is  accurate,  we  have  yet 
to  verify  one  crucial  feature  of  that  diagram.  Specifically,  the  argument  assumed  that  if  we 
enter  region  II  at  the  point  P  and  travel  within  this  region  keeping  v'  >  0  and  u'  <  0,  then 
we  cannot  reach  region  III.  (Clearly,  no  problem  arises  from  the  possibility  of  re-entering 
region  I.).  Thus  we  have  to  show  that  the  situation  shown  in  Diagram  2  is  impossible. 
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Diagram  2 
Not  possible 
No  such  Q  exists 
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Analytically,  we  must  show  that 


^  ^+('>^2)  for  u,  <  ui  <  112  Jeq7io(3.16) 

Defining  Hi  —  —H(ui)  and  Ai  =  —Am^H"(ui)H{ui),  then  (3.16)  boils  down  to 

(3.17)  -2/7i{l  +  v^l  -  Ai}Mi  <  -27/2(1  -  ^J^^A:)|A2  . 

Now,  the  convexity  of  /  (and  hence  of  H)  together  with  the  fact  that  H  vanishes  at  u±  imply 
that  Hi  >  H2  >  0;  these  facts  together  with  our  assumption  (3.3)  imply  that 

(3.18)  1  >  A-  >  0  . 


Therefore,  a  sufficient  condition  for  (3.17)  to  hold  is  that 

(3.19)  {1  +  —  Ai}/Ai  >  {1  —  \Jl  —  A2)I A2 


for  all  Ai  satisfying  (3.18).  A  little  algebraic  manipulation  shows  that  this  is  indeed  the  case. 
Consequently,  (3.16)  holds,  i.e.,  no  point  such  as  the  point  Q  in  Diagram  2  can  exist,  and  so 
the  argument  based  on  Diagram  1  is  valid. 

As  m  increases  past  the  value  that  makes  equality  hold  in  (3.3),  we  obtain  the  situation 
shown  in  Diagram  3. 

Namely,  a  gap  appears  in  region  //,  through  which  our  trajectory  might  possibly  plunge 
into  the  abyss  of  region  III.  Hence  we  cannot  say  whether  a  trajectory  joining  u~  to  u"*" 
exists  or  not.  Finally,  when  m  increases  past  the  value  that  makes  equality  hold  in  (3.4), 
then  the  phase-plane  looks  like  Diagram  4,  and  the  descent  of  our  trajectory  to  —00  becomes 
a  certainty.  □ 

We  close  this  section  by  quantifying  Rosenau’s  statement  [R]  that  the  travelling-wave 
solutions  of  the  R-C-E  equation  (1.1)  give  narrower  shock  layers  than  those  of  the  viscosity 
equation  (1.4).  We  will  adopt  as  our  measure  of  shock  width,  w  =  (u_  —  u+)/|u',|  with  u', 
evaluated  at  the  point  ii»  at  which  H'{u)  =  0.  (It  should  be  noted,  however,  that  although 
this  value  of  u'^  is  always  meucimum  for  (1.4),  it  will  be  maximum  for  (1.1)  only  if  the 
trajectory  does  not  enter  region  II  of  Diagram  1,  which  in  turn  is  guaranteed  by  our  analysis 
only  when  the  curve  v  =  v+{u)  has  its  unique  local  minimum  at  u,.)  Since  the  relevant  value 
of  u'^  for  (1.4)  is  given  by  u'^  =  —H[u,),  while  the  value  of  u\  for  (1.1)  is  u.^(u,),  the  estimate 


(3.20) 


"^Chapman — Eiwkog 
^viscou* 


<  1 


follows  from  the  simple  lemma:  If  a  quaaratic  equation  has  real  roots,  then  the  root  closer 
to  zero  lies  between  r  and  2r,  where  r  is  the  root  of  its  linear  part. 
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4.  Entropy  Solutions  and  the  Zero  Mean-Free-Path  Limit 

The  parameter  m  does  not  play  a  role  in  our  analysis  in  this  section,  and  so  will  be  set 
equal  to  1  for  convenience. 

Since  solutions  of  the  R-C-E  equation  (1.1)  may  contain  singularities  weak  solutions  must 
be  admitted.  Since  the  latter  need  not  be  unique,  we  single  out  an  “entropy”  solution  of  the 
R-C-E  equation  (1.1)  as  the  one  satisfying  the  Krushkov-like  [K]  inequality 

dt\u,-c\  +  5*[spn(Ue  -  c){f{Ue)  -  /(c)}] 

(4.1) 

£  -;{!“■  “  “^1  ~  >«'*(“•  -  '=)Q-  •  (“■  - 1:)}  . 

for  all  real  c’s.  In  particular,  by  choosing  c  =  -f  sup  [u*]  (respectively,  c  =  —  sup  |u*|),  we 
obtain  from  (4.1)  that  u*  is  a  supersolution  (respectively,  a  subsolution)  of  (1.1),  and  hence 

(1.1)  is  satisfied  in  the  sense  of  distributions.  We  turn  to  show  that  (4.1)  admits  a  unique 
solution  Ug,  and  that  this  solution  converges  to  the  u  *'  ue  entropy  solution  of  (1.4)  as  e  goes 
to  zero. 

Theorem  4.1.  For  any  uq  in  BV  there  exists  a  unique  solution  Ug  of  the  R-C-E  equa¬ 
tion  (4-1))  as  £  i  0,  Ug  converges  in  to  the  unique  entropy  solution  of  the  inviscid 
conservation  law  (l.f)- 

Proof. 

Add  the  artificial  viscosity  term  to  the  right  side  of  (1.1);  the  resulting  equation  has 
a  unique  smooth  solution  u\.  By  a  straightforward  adaptation  of  Krushkov’s  [K,  section  4] 
proof  for  the  artificial  viscosity  method  for  (1-4),  we  obtain  that  the  set  {Ug}f>o  is  bounded 
in  BV  (uniformly  in  e  and  6)  and  precompact  in  L^,  and  hence  that  a  subsequence  converges 
as  5  —>  0  to  a  solution  Ug  of  (4.1).  Similarly,  by  the  argument  on  pages  224-5  of  [K]  we 
obtain  from  (4.1)  the  consequence 


noo 

{[Ue  -  -t-  3gn{ug  -  Vg)[f{ug)  -  /(u,)]$x}da 

'OO 


(4.2) 


dt 


noo  1 

-{[tie  -  Vg\  -  sgn{ug  -  Vg)Qg  *  (uj  -  Vg)}^dxdt 

-oo  £ 


where  $  is  an  arbitrary  nonnegative  test  function. 

Next,  we  remark  that  the  expression  inside  the  curly  brackets  on  the  right  of  (4.2)  need 
not  be  positive,  but  in  view  of  (2.5),  its  spatial  integral  is.  Therefore,  by  choosing  $(t,x)  = 
4^i(<)$2(x),  and  letting  $2  l®nd  to  the  function  that  is  identically  one,  we  obtain 

(4.3)  [  f  |u,  -  >  0  . 

Jo  J —00 
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Continuing  as  in  [K],  we  let  approach  the  indicator  function  of  the  interval  [0,t]  to 
conclude 


(4.4)  \\u,{t)  -  u.COIIli  <  \\Uc{0)  -  Ue(0)llii  . 

In  particular,  this  shows  that  the  solution  of  (4.1)  is  unique. 

The  solutions  {u*}  of  (4.1)  inherit  the  BV  bound  of  the  {u,},  and  the  argument  of  section 
4  of  [K]  shows  that  this  bound  implies  precompactness  in  L^.  Hence  as  e  — »  0  a  subsequence 
converges  to  a  weak  solution  u  of  (1.4).  Because  the  right  side  of  (4.2)  is  known  to  be  positive 
only  when  $  has  no  dependence  on  x,  we  cannot  use  the  entropies  of  (1.4)  (as  in  [Ta])  to 
conclude  that  u  is  the  entropy  solution  of  (1.4).  However,  (4.4)  implies  the  corresponding 
estimate  for  the  weak  solutions  u  and  v  obtained  in  the  limit  as  e  goes  to  zero,  and  by  an 
argument  of  Lax  [L]  this  suffices  to  show  that  we  obtain  the  entropy  solution:  It  is  not  hard 
to  see  that  when  (1.4)  has  a  smooth  solution  then  our  scheme  must  converge  to  that  solution. 
Hence  by  the  corollary  to  theorem  (3.5)  of  [L],  any  solution  u  of  (1.4)  obtained  in  the  limit 
e  — >  0  from  (4.1)  has  the  property  that  all  of  its  discontinuities  satisfy  the  generalized  Lax 
shock  inequalities.  By  theorem  (3.5)  of  [L],  this  implies  that  u  is  the  unique  entropy  solution 
of  (1.4).  Finally,  since  any  sequence  of  es  tending  to  zero  has  a  convergent  subsequence,  the 
uniqueness  of  the  limit  shows  that  convergence  holds  without  passing  to  a  sequence.  □ 


5.  The  Convergence  Rate  of  the  Zero  Mean-Free-Path  Limit 

Theorem  4.1  shows  that  the  R-C-E  equation  (1.1)  retains  several  properties  of  the  viscous 
conservation  law  (1.2).  In  particular,  (4.4)  asserts  that  the  solution  operator  is  an  L^- 
contraction,  and  hence  by  conservation  plus  translation  invariance  it  is  monotone  [CM, 
Lemma  3.2],  and  by  translation  invariance  it  is  BV-bounded: 

(5.1)  ||ue(0llBV  <  |K(0)||bv 


Next  we  show  that  the  nonlinear  R-C-E  equation  (1.1)  also  satisfies  Oleinik’s  E-entropy 
condition,  e.g.,  [Sm],[T]. 

THEOREM  5.1.  Assume  f"'>a>  0.  Then  the  following  a  priori  estimate  hold^ 

1 


(5.2) 


'^e(0llLtp+  — 


ll^*(0)llz,i+  + 


,  t  >  0. 


*Wc  let  ll^llLtpi  ||i^||i,<p+  and  denote  respectively  eB85up*^y 

and 


x-v 


esssupx^y 


i-V 


14 


Remark.  The  inequality  (5.2)  implies  that  the  positive-variation  and  hence  the  total- 
variation  of  Ue{t)  decays  in  time.  Furthermore,  this  gives  us  another  proof  of  the  zero 
mean-free-path  convergence  to  the  entropy  solution  of  (1.4)  for  any  I-j^-initial  data  uq  (cf. 
Corollary  5.2). 

Proof.  We  add  the  artificial  viscosity  term  6uxx  to  regularize  (1.1),  obtaining 

(5.3)  a,ni  +  =  -4-{uf  -  (?„,  .  ui)  +  Wjuf. 

Differentiation  of  (5.3)  yields  for  w  =  dxU^, 

(5.4)  dtw  -I-  f‘(ul)d:,w  -f  /"(uf)u;^  == - *  u;}  -f  6dlw. 

Hence,  since  /  >  a  >  0,  it  follows  that  W(t)  =  maXxW^t)  is  governed  by  the  differential 
inequahty 

(5.5)  W(t)  -f  aW\t)  <  *  W}  <  0 

and  (5.2)  follows  by  letting  S  I  O.Q 

Theorem  5.1  shows  that  solutions  of  the  R-C-E  equation  (i.i)  are  Zzp '^-stable.  Moreover, 
(5.1)  implies  that  the  Zzp'-size  of  their  truncation  if  of  order  0{c),  for 

(5.6)  \\dtu,  +dxf{u,)\\up'  =  e\\Qmc  ll«e(0llBV  <  el|ue(0)||Bv. 

Using  the  result  of  [T]  we  conclude  that  the  Lip'-  convergence  rate  of  the  R-C-E  solutions 
to  the  corresponding  entropy  solution  is  also  of  order  0{e). 

Corollary  5.2.  Assume  that  f  >  a  >  0,  and  let  Ug  be  the  unique  R-C-E  solution  of  (^.1) 
subject  to  initial  conditions  Ug^O)  =  ^  thenUg  converges  to  the  unique  entropy  solution 

of  (1.4)  OLTid  the  following  error  estimates  hold 

\\ug{t)  —  u{t)\\w-‘.r  <  Const  ■  ,  l<p<oo,  5  =  0,1. 

Remark.  The  choice  (s,p)  =  (1, 1)  corresponds  to  Lip'-convergence  rate  of  order  0(e).  The 
choice  (s,p)  =  (0, 1)  corresponds  to  the  usual  viscous  LUconvergence  rate  of  order  O(e^). 
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which  motivate  the  present  study. 
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